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We describe a recent implementation of the combined GW and dynamical mean field (DMFT) 
method "GW+DMFT" for the two-dimensional Hubbard model with on-site and nearest-neighbor 
repulsion. We clarify the relation of the GW+DMFT scheme to alternative approaches in the liter- 
ature, and discuss the corresponding approximations to the free energy functional of the model. 
Furthermore, we describe a numerically exact technique for the self-consistent solution of the 
GW-I-DMFT equations, namely the hybridization expansion continuous-time algorithm for solv- 
ing the dynamical impurity model that arises within the GW-I-DMFT scheme. We compute the 
low-temperature phase diagram of the extended Hubbard model, addressing the metal-insulator 
transition at small intersite interactions and the transition to a charge-ordered state for stronger 
intersite repulsions. Within the GW-I-DMFT framework, as in extended DMFT, the intersite re- 
pulsion translates into a frequency-dependence of the local effective interaction. We analyze this 
dependence by extracting a charateristic plasma frequency and show how it affects the local spectral 
function. 

PACS numbers: 71.10.Fd 



I. INTRODUCTION 



Understanding the effects of strong electronic correla- 
tions in lattice systems remains a challenge in condensed 
matter physics. The competition between the tendency 
of electrons to delocalize due to the corresponding gain 
in kinetic energy and localisation by the Coulomb in- 
teraction gives rise to a panoply of interesting phenom- 
ena, ranging from simple mass enhancements in the sense 
of Landau theory to charge-, spin- or orbital ordering 
phenomena. To reduce the complexity of the problem, 
while still keeping the main qualitative effects, e. g., of 
the delocalisation-localisation transition, one often re- 
sorts to low-energy effective models such as the Hub- 
bard or Anderson lattice models. The two-dimensional 
single-band Hubbard Hamiltonian with a static on-site 
repulsive interaction U for example is believed to cap- 
ture the physics of the high-tcmpcrature superconducting 
cuprates»i Charge-ordering transitions can be captured 
when also an intersite interaction term, mimicking the 
longer range Coulomb interactions, is retained. This mo- 
tivates the construction of the extended Hubbard model, 
where charge-ordering effects and screening of the local 
interactions due to the non-local ones are included in ad- 
dition to the pure Hubbard model physics. 

In the paramagnetic phase at half-filling, the Hubbard 
model exhibits a Mott transition from a metal to a Mott 
insulator whose spectral function is characterized by a 
gap at the Fermi energy and Hubbard bands correspond- 
ing to atomic-like excitations. This behavior is captured 
by computational schemes such as the Dynamical Mean 
Field Theory (DMFT) (see Refs. and i for review) or 



its extensions (C-DMFT^i^, DCA^, dual fermions^, dual 
bosons*^, DTA^, BMFT+T,!^). A formalism which al- 
lows to treat screening by non-local interactions is ex- 
tended DMFT (EDMFT) Jii-i^ Its combination with the 
ab-initio GW approacbi^^— introduces some momentum- 
dependence also into the self-energy, capturing the inter- 
play of screening and non-local correlations. This scheme 
allows for a self-consistent computation of the effective 
"Hubbard U" in a solid and a fully parameter-free ab- 
initio simulation approach for correlated materials. The 
idea is to take the local part of the self-energy from the 
EDMFT calculation and supplement it by the nonlocal 
component of the GW self-energy. A rigorous functional 
formulation - which is detailed in section II - puts this 
theory on the same level of mathematical rigor as, e. g., 
the functional formulation of Hohenberg-Kohn density 
functional theory. 

Despite these promising perspectives, the technical dif- 
ficulties associated with the numerical treatment of the 
frequency-dependent effective interaction have so far pre- 
vented a self-consistent calculation of spectral properties 
within GW-I-DMFT, even at the model level. Recent 
progress, both within approximate schemes^ and nu- 
merically exact Monte Carlo techniques?^ is currently 
giving a new impact to the field i^"— In particular, the 
development of efficient continuous-time Monte Carlo 
techniques ^i^^ generalized to dynamical interactions in 
Ref. [20, now allows for the fully self-consistent solution of 
the GW-I-DMFT equations, with high enough accuracy 
to also extract spectral properties. 

In this paper we use these state-of-the-art techniques 
to study an extended Hubbard model with on-site inter- 
action U and nearest-neighbor repulsive interaction V, 
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and explore the interplay of screening and non-local cor- 
relations. We show that this model, if solved within the 
EDMFT framework, captures dynamical screening effects 
related to the nonlocal interaction V, high-energy satel- 
lite features in the one-particle spectra and, for large V, 
a transition to a charge-ordered state. We then proceed 
to study how diagrammatic corrections to EDMFT in 
the form of momentum-dependent GW contributions to 
the self-energy modify this picture, and compare the re- 
sults of self-consistent and non-self-consistent implemen- 
tations. 

The paper is organized as follows: in Section II, we 
discuss the model. In particular, we will show how the 
Hamiltonian formulation of the problem can be recast 
into an action or functional formulation. We will ex- 
plicitly construct two flavors of free energy functionals 
whose stationary points would give the exact solution of 
the model, the first one due to Almbladh et al^ (see also 
Ref. 113), the second one constructed by Sun et al.i^ In 
section III, we present different methods of solution: ex- 
tended DMFT, GW and combined schemes. We argue 
that the GW-FDMFT formalism can be derived from the 
Almbladh free energy functional, and discuss the differ- 
ences to the scheme proposed in Ref. [13. The latter 
stems from another energy functional, which would cor- 
respond to a "GD-^SOPT-I-DMFT" method, where D is 
the boson propagator associated to screening of the non- 
local interaction only, used in a GW-like fashion, sup- 
plemented by second order perturbation theory (SOFT) 
for the non-local effects of the local interaction and dy- 
namical mean field theory. Numerical techniques for the 
solution of the equations are described in Section IV and 
the computational scheme used in the present work is 
summarized in Section V. The results of our study are 
presented in Section VI. Section VII contains a summary 
of our most important findings, and provides perspectives 
how our work inserts itself into the field. 



II. MODEL 
A. Hamiltonian Formulation 

We consider the single-band U-V Hubbard model on 
a two-dimensional square lattice, defined by the grand- 
canonical Hamiltonian 

H = -tY^ {c\^cj„ + h.(^ - ^x'Yn, 

{tj)a i 

+ U'^Hi^riii + vYninj, (1) 

where Ci^ and denote the annihilation and creation 
operators of a particle of spin a ~tii at the lattice site 
i, ni„ = c\^Ci^, and Ui = + n^. denotes the 

sum over nearest-neighbor bonds, t > is the hopping 
integral between two neighboring sites, /i is the chemical 



potential, U the on-site interaction between electrons of 
opposite spin and V the interaction between two elec- 
trons on neighboring sites, irrespective of their spin. The 
number of nearest neighbors is z = 2(i = 4, where d is 
the dimension. 

This model can be obtained from first principles in 
realistic calculations, as discussed in Appendix [X] The 
limiting case = corresponds to the conventional Hub- 
bard modelj^"— In the present study, we will limit our- 
selves to repulsive interactions, U > Q and V > Q. 

In the limit of large the extended Hubbard model 
has been shown to display a transition to a charge- 
ordered state characterized by a freezing of charge car- 
riers and a spatial modulation of the charge density 
(charge-density wave, CDW) . This can be explained by a 
simple energetical argument in the strongly correlated 
regime [U ^ t): while for U /V ^ 1, electrons will 
lower their energy by arranging themselves one per site 
to minimize the on-site repulsion, for V/U ^ 1, elec- 
trons will minimize their off-site repulsion by choosing 
an arrangement such that one sublattice is occupied by 
two electrons per site while the other is empty, leading 
to a commensurate charge order. In the metallic phase 
([/ ^ i), the effect of V is more easily understood in 
terms of screening: the charge fluctuations induced by 
the V term lead to a reduction of the local effective in- 
teraction. 

The U-V Hubbard model has been studied by a va- 
riety of approximate methods. An early study in the 
zero-overlap limit {U /t ^ 1) has predicted a phase tran- 
sition between a Mott insulator and a charge-ordered 
(CO) insulator at = J7/z at zero temperature,'^^ while 
weak-coupling mean-field studies have predicted a tran- 
sition between antiferromagnetic order (AFM, i.e. com- 
mensurate SDW) and charge order {i.e. commensurate 
CDW) at the same boundaryj^i^ This has been con- 
firmed by Monte-Carlo calculations in two dimensions!^ 
The Vc = U / z boundary has been shown to hold in the 
U /V — limit as well as in the U /V — )■ oo limit by 
a study at half-filling in the infinite dimensional limit<^ 
Higher-order corrections have been considered in Refs.[36l 
and [37I leading to a renormalization of the critical tem- 
perature and order parameter, as well as the predic- 
tion of phase separation in the zero-temperature limit. 
More recently, variational cluster— and two-particle self- 
consistent approaches'^^ have been applied to the U-V 
Hubbard model. 

A first DMFT treatment described the opening of a 
"pseudo-gap" and the reentrant behavior of the critical 
Vc as a function of temperature.'**' In this simple scheme 
the V-ierm only contributed at the Hartree level by shift- 
ing the chemical potential, since in the limit of infinite 
dimensions, the contributions beyond Hartree of non- 
local terms vanish, while local terms such as the on-site 
Hubbard U do not vanish under rescalingi^ The screen- 
ing effects contained in the T^-term are not captured 
by standard single-site DMFT. Therefore, an extended 
DMFT (EDMFT) scheme has been proposed to remedy 
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this shortcoming. -'^^"^^ Within this scheme, the non-local 
interactions induce a frequency-dependence of the effec- 
tive local interaction and lead to a sizeable reduction of 
the static value of U. In addition, Refs. TJ and^ showed 
that one of the consequences of adding a spatially non- 
local contribution to the self-energy is to make the system 
more insulating. However, none of the previous studies 
has given a precise account of how U and V affect the 
properties of the local frequency-dependent interactions 
nor how the latter in turn modify the spectral properties 
of the system. A systematic investigation of the effect 
of non-local GW contributions to the self-energy is also 
lacking. 

In the present work, we try to clarify these issues, re- 
stricting ourselves to the paramagnetic phase at half- 
filling. We will give all energies in units of the half- 
bandwidth, D = At. The inverse temperature will be 
denoted by /3 = 1/T {kg = 1). 



B. Action Formulation 

The solution of model ([T]) amounts to computing the 
Green functions and other correlation functions. For this 
purpose, it is convenient to write the grand-canonical 
partition function Z — Tre~^^ as a coherent-state path 
integral^ Z = J 'D[c*,Ci]e~^ where 

S[c*,c] = J ^ cKt) iidr - m) S., + t,,) c,„{t) 

+UY,n,^{T)n,^{T) + i 5]<n,(r)n,(r) ) . (2) 



c* and Ci denote conjugate anticommuting Grassmann 
fields for site i {ci{T + (3) = ~Ci{T)), tij — —td^ij), 
v^j = V^^(ij) and S(^ij-^ = 1 if i and j are nearest neighbors 
and otherwise, r is the imaginary time variable. We 
will denote by a;„ = (2n -I- 1)tt/ (3 (y^ = 2nn/ (3) the cor- 
responding fermionic (bosonic) Matsubara frequencies. 
The Fourier transforms of tij and vf^ on the lattice are: 



jjL + U /2. The action hence becomes 



-~2t (cos{knc) + cos{ky)) , 
2V {cos{kx) + cos(fcy)) . 



(3) 
(4) 



Using the identity niTii — 2ni^ni^ + JT-i, we can rewrite 
the interaction terms of Eq. ^ as \ '^ij'Vijni{T)nj{T) 
with = U Sij + VS(^ij^ , or 



Vk = U + vf, (5) 



provided that we shift the chemical potential fi ^ fl 



S[c,c] ^ J dr^Yl c«<t(t) {{dr - A) % + t,j) Cj„{t 



C. Functional Formulation 

The problem of finding the solution to the Hamilto- 
nian model ([Ij or calculating the Green's function cor- 
responding to the action ([2|) can equivalently be formu- 
lated in a functional language. The familiar Luttinger- 
Ward or Baym-Kadanoff functionals provide examples of 
such a construction. In the present context, a formula- 
tion in terms of the free energy written as a functional of 
both, the Green's function G and the screened Coulomb 
interaction W, is the method of choice, since the com- 
bined GW-I-DMFT method can naturally be viewed as a 
specific approximation to such a functional. Indeed, the 
GW+DMFT solution as formulated in Refs. [H and [13 
can be derived as a stationary point (G, W) of the free 
energy functional introduced by Almbladh et al.r^ after 
approximating the correlation part of this functional by a 
combination of local and non-local terms stemming from 
DMFT and GW, respectively. To draw an illustrative 
analogy, GW-I-DMFT provides an approximation to the 
correlation part of the Almbladh free energy functional, 
on the same footing as the local density approximation 
of density functional theory^ is an approximation to the 
exchange-correlation part of the Hohenberg-Kohn func- 
tional of the energy. 

In the literature, several variants of functionals of G 
and W have been discussed r^iiii^ and different deriva- 
tions given. Using a Hubbard-Stratonovich (HS) decou- 
pling as in Ref. Ilg . we discuss two fiavors of free en- 
ergy functionals which differ by the choice of the part of 
the interaction that is decoupled by the HS transforma- 
tion. The first one reproduces the 'J functional intro- 
duced by Almbladh on which the GW-I-DMFT construc- 
tion of Ref. 116| is based. The second one is a variant that 
was used in the study of the extended Hubbard model in 
Ref. M- 

The Hubbard-Stratonovich transformation!^ relies on 
the following identity: 



exp 



■D[xi{t),X2{t), ...] 
^(27r)^detA 



xexp^y" dr i^-^x^{T)[A %Xj{t) T Xi{T)b,{T)yj{7) 

where A is a real symmetric positive-definite matrix, 
6i(r) and Xi{T) are periodic real fields (5i(r-|-/3) = ?>i(r)), 
and summation over repeated indices is understood. In 
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the foUowing, we will choose the upper sign for the last 
term. 

Decoupling the whole (local and non-local) interac- 
tion ^'^ij Vijni{T)nj{T) by the HS transformation cor- 
responds to applying the above formula ([7]) to the inter- 
action term in Eq. ([6]), that is to the choice bi = iui, 
Aij = Vij and Xi = ipi. This choice - denoted as HS- 
UV in the following - leads to the construction of the 
^ functional as in Almbladh et al.^^ and the formal- 
ism of Ref. [TtI . Within the approximation introduced 
in the next section this corresponds to the GW-(-DMFT 
scheme as introduced in Ref. Ilg . This approach (which 
was also used in Ref. flsh relies on the argument that 
the two terms, which represent different matrix elements 
of the same interaction, should be treated on the same 
footing. 

In Ref. [13, on the other hand, the HS transform 
has been applied only to the non-local interaction term 
^ t;g'ni(T)nj(T) in the action ([2]). This approach, 
dubbed HS-V in the following, leads to a modified free 
energy functional, which we denote by ^'y. It was moti- 
vated by the aforementioned fact that in the limit of in- 
finite dimensions, with the non-local interaction rescaled 
as V ^ V/z, the non-local term results in a trivial shift 
of the chemical potential, while the on-site interaction 
remains non-trivial, justifying a separate treatment for 
the non-local term. 

We will first explicitly work out the construction of 
the two functionals ^E* and ^I'v and in section IVIl we 
will compare the results from both decoupling strate- 
gies. The results in section IVI Bl finally, are based on 
the HS-UV decoupling, that is, on the ^ functional and 
the GW-hDMFT formalism of Ref. M- 



1. "UV-Decoupling" : the ^' Functional 

In the HS-UV decoupling scheme, the full interac- 
tion term is decoupled via an auxiliary bosonic field (pi. 



This has an important consequence: even a first-order 
diagram in this new interaction contains diagrams of in- 
finite order in the electron-electron interaction. Making 
even simple approximations on the new action can thus 
lead to non-trivial diagrams for the original action. More- 
over, the electron-boson vertex icpiUi is local. This local- 
ity ensures that in the limit of infinite dimensions, the 
interactions (and in particular V) will contribute beyond 
the Hartree level. 

The generating functional of correlation functions is 
obtained by introducing the bilinear sources Jf and Jb, 
coupling to the fermionic and bosonic operators respec- 
tively, so that the action becomes S[c*,c,4>] — S[Ji,Jh], 
with 

5[Jf, Jb] =Y, dTdT'(^Jf^,,iT,T')c*iT)c,{T') 



+ ^Jb.vir,T')Mr)<P,{r')]. (9) 



The fermionic and bosonic Green's functions for this ac- 
tion are G^,{t - r') = -^{Tc^{T)c]iT')) = Sn/SJujir, r') 
and W,,iT-T') = {T(j),{r)^j{r')) = -2Sn/SA,,j{T,T'), 
where we have defined 



lnZ[Jf, Jb] = -lnTre-^["*'"''^l+^l-^f'^-l. (10) 



The non-interacting Green's functions (obtained by 
setting icpiUi — 0) are respectively G\a=Q = Gq and 

W{k,ii/n)\a=o = [f^k^] = ^fc- This gives a first hint as 
to the physical meaning of W: without renormalization 
by the auxiliary bosons, it corresponds to the bare inter- 
action. Coupling to the bosons, which represent density 
fiuctuations of the system, explicitly introduces screening 
into the physical description. 

We next perform a Legendre transformation with re- 
spect to the sources Jf and Jb, 



Choosing bi = iui, Ai 



and Xi 



i, the transfor- 



mation (O applied to the action ([2]) leads to 



T[G,W]^n[Ji[G],MW]] 



Tr Jf G 4- ^Tr JbVF (11) 
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c ,c, ( 



1 



ij i ) 

(8) 

where we introduced the fermionic lattice Hartree 
Green's function G^ defined by [G^^^],^- = {{-dr + 

IJ. + For later use, we have moreover inserted 

a coupling constant a in front of the fermion-boson cou- 
pling term. The physically relevant case corresponds to 
a=l. 

The HS transformation replaces an electron-electron 
interaction by an electron-boson interaction, and intro- 
duces a new variable, the auxiliary real boson field <j). 



E 

.-ll 



with the reciprocity relations Ji = and Jb = 2^f. 
The physical Green's functions will be obtained by set- 
ting Jf = and Jb = 0, or equivalently, by requiring the 
stationarity of F with respect to G and W . 
The free energy functional F can be written as 



Fq=i — Fq=o 



where we have defined 



Jo 



da 



(12) 



(13) 



F is the well-known Baym-Kadanoff functional,''^ while 
is the extension of the Luttinger-Ward functional $[G] 
to one- and two-particle propagators^^ 
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The non-interacting {a = 0) part of the F-functional 
is readily evaluated as 

ra=o = Trln(-G') - Tt{G^^ - G-^)G 

-^TrlnW +^Tt{v-^ -W'^)W. (14) 

Indeed, when a = 0, the action becomes Gaus- 
sian and thus explicitly integrable, namely: ^a=o — 
-lnDet[-Go ^ + Jt] - ln(Det[i;-i - Jb])^/^ The above 
definition G = Sil/SJf imposes (Go ^ - Jf)G = 1 and 
similarly {v^^ — Jh)W — 1 yielding Eq. ([I4|) . Finally, 

SG 



stationarity of the full F reads ^ — — 



SG 

^ for G and = -^{W-^ - y-^) + |^ 



^ SG 

for W. Defining the self energies as 



(5* 

SG' 



n 



yields Dyson's equations for G and W: 



G — Gfi 



n. 



(15) 



(16) 



Being "^'-derivable" , these self-energies will obey conser- 
vation rules. 

The above formulation shows that, formally, solving 
the lattice problem defined by Eq. ([T71) amounts to eval- 
uating the corresponding 'i'-functional, from which S and 
n, and in turn G and W can be derived. In section Hill 
we will describe two complementary ways of approximat- 
ing this functional, EDMFT and GW, before showing 
how to merge the two approaches, thus arriving at the 
GW+DMFT free energy functional. 



2. "V-Decoupling" : The 'i'v Functional 

In the HS-V scheme, proposed in Ref. [13, only the 
non-local interaction term is decoupled via an auxiliary 
bosonic field Choosing bi = im, Aij = v^j and Xi = 

, the transformation JT]) applied to the action ^ leads 

to 



S 



0] = f drl - ^ cUt) [Go 1] c,,{t) 

I ija 

+ aU^ mi (T)ni^ (r) I 

(17) 

where we introduced the non-interacting fermionic lattice 
Green function Go defined by [Gq ^] ■ - = ((— 3r + m)% — 
tij). Again, a coupling constant a was introduced, and 
the physical case corresponds to a=l. Now, however, the 
coupling constant is not only a switch for turning on or 



off the fermion-boson coupling but at the same time also 
the local Hubbard interaction. 

In principle, the interaction should be a positive def- 
inite matrix in order for the Gaussian integrals invoked 
in the HS transformation to converge. Unlike the situa- 
tion in the HS-UV decoupling where U and V are matrix 
elements of the screened Coulomb interaction, which is 
positive definite, this is not the case for the interaction 
of HS-V, v^j. This issue can be dealt with by adding 
an auxiliary identity matrix multiplied by a large enough 
constantJ^ In practice, however, the simulation results 
are not affected by the value of this constant. 

As before, the generating functional of correlation 
functions is obtained by introducing source terms. The 
fermionic Green's functions for this action is unchanged 
compared to the UV-decoupling case: Gij{T — r') — 
~{Tci{T)cj{T')) = 6n/6Ji^ij{T,T'). The bosonic propa- 
gator formally still reads (r — r') = {T(j)i{T)(j)j{T')) = 
—26^}/SJh,ij{T,T'). It does not, however, correspond to 
the screened interaction, as in the HS-UV scheme: in 
the case of vanishing fermion-boson coupling, the bosonic 
propagator reduces by construction to only the non-local 
part of the bare interaction. 

The construction of the free energy functional F pro- 
ceeds as before by Legendre transformation with respect 
to the sources Jf and Jb, 

Tv[G,D] = n[Jt[GlMD]] - TrJfG+ ^TrJb^, (18) 

with the reciprocity relations Jf = — and Jb — S^Jj^. 
The physical Green's functions will be obtained by set- 
ting Jf — and Jb — 0, or equivalently, by requir- 
ing the stationarity of Fy with respect to G and D. 
Thanks to the choice of the coupling constant a in front 
of the interaction and boson-fermion coupling terms, 
a (U '^ni-^rii^ + i'^^(j)ini), Fy acquires the same form 

as before, Tv,a=i = ^v.a=o + ^'v, with = da^, 
but it is now a functional of G and D. 

The non-interacting {a — 0) part of the F-functional 
reads 

ry,a=o = Trln(-G) - 1^^{G^^ - G-^)G 

-^Trlniy + iTr((w"')-i - D-^)D. (19) 

Finally, stationarity of the full Ty reproduces the 
fermionic Dyson equation for the Green's function and 
self-energy. For the bosonic part, however, we obtain 
= -i(L»-i - (w"')~^) + ^ for £1. The bosonic self- 
energy 



Hi, 



' SD ' 



(20) 



is thus not equal to the physical polarization of the sys- 
tem. 

Again, solving the lattice problem defined by Eq. p7)) 
amounts to evaluating the corresponding ^'y-functional. 
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from which S and Ily, and in turn G and D can be de- 
rived. Compared to the previous case of the ^'-functional, 
however, the sublety of D being the screened non-local 
interaction (not equal to the full W) requires additional 
care in the construction of a combined DMFT scheme. 



III. METHODS OF SOLUTION 

A. EDMFT 

^E" is a functional of the fermionic and bosonic 
Green functions Gy and Wij. EDMFT approximates 



This action is very similar to the lattice action (Eq. ()17|) '). 
with Go replaced by an appropriately defined dynamical 
Q decribing the excursions of an electron in the lattice 
from a given site (the impurity, defined by operators 
and c^) and back, and the bare and instantaneous in- 
teraction V (or u"' in HS-V) replaced by the retarded 
interaction U (or 2?) . The effective actions (|2T|) and (|22|) 



-lEDMFT _ 
'eff, HS-UV — 

qEDMFT _ 
'-'eff, HS-V — 



which feature a retarded retarded interaction U{t — t') 
(for HS-UV) or V(t - t') (for HS-V) between charges. 

The solution of this impurity problem, described in 
detail in Section |TVl requires the calculation of the one- 
particle Green's functions Gioc = —{Tc{t)c^(Qi)) and 
Wioc = {T4>{t)4>{Q)) (or Aoc for HS-V). From Gioc 
and W^ioc, one computes the corresponding self-energies, 
Sloe = Q-^ - Gi-J and Hioc = U-^ - W{-^1 (or (Hv)ioc = 



^[G,j,Wrj\ by *^^*'^^^[Gii,Wii] = ^[G,A,,W,A,]. 
Similarly, [G,„ Ai] = ^v[Gh5,j,D,Aj\. The 

local Green functions can be obtained by solving an aux- 
iliary effective local problem defined by the action 



(21) 



(22) 



are obtained by integrating out all sites but one in the 
lattice action and taking the infinite-dimensional limit. 
The derivation of the action for the HS-UV scheme, as 
well as the EDMFT loop sketched below, are presented 
in Appendix [BJ 

Integrating out the 0-field in Eqs. (PT|) and yields 
the impurity actions 



(23) 

2 

(24) 



- D^^l). The EDMFT approximation consists in 
identifying the impurity self-energies with the lattice self- 
energies: S;(fc, iwn) Eioc(ia;„), H(fc,w„) w Il\oc{iVn)- 
This allows one to evaluate (approximate) lattice Green 
functions G(fc, iw„) and W{k, ivn) (or D{k, ivn)) through 
Dyson's equation and estimates of the local lattice Green 
functions by summation over k. Eventually, one obtains 
updated Q and U (or V) via 



JO 



dTdT'(l){T)U-^{T-T')cl){T)+i / dT(l)(T)n{T) 

Jo 

13 r0 



S^s%V = - / dTdr'Y,ci{T)g"\T-r')c,{T')+ f drUn^{T)n,{T) 
Jo „ Jo 



dTdT'(l){T)V-^{T -T')(l){T')+i I dT<j>{T)n{T). 

Jo 

I 



nP T nP 1 

/ dTdT'y24iT)S^'^i'r~T')(^'r{T') + - dr dr' n{T)U {t ~ t' )7i{t' ) - -Tv \nU , 
Jo „ 2 Jq 



2 

P 



rP rP T rP 

■ / drdr' V cl{T)g-^{T - t')c„(t') + / dTUn-^{T)ni{T) + - dTdT'n{T)V{T ~ r')n(T') 
Jo „ Jo 2 Jo 



2?-i = 



A"oc[(ny)ioc 



5j1oC7 

IIioc , 

] + (ny)loc. 



(25) 

(26) 
(27) 



B. The GW approximation 



While EDMFT can treat strong local correlations, it 
completely neglects the nonlocal contributions to the self- 
energy. A complementary approach, which treats spatial 
fluctuations, but works reliably only in the weakly corre- 
lated regime is the GW method.""* The GW approxima- 
tion has been used extensively to investigate the prop- 
erties of weakly correlated materials, such as the band 
gaps of semi-conductors (see e.g. Ref . l49l) . In these ma- 
terials GW correctly accounts for the screening effects of 
the electrons at the random-phase approximation (RPA) 
level. Schematically, for a general Coulomb interaction 
v{r) ~ 1/r, one-shot GW consists in replacing the bare 
interaction v of the Fock self-energy Yip ^ Gqv by the 
screened interaction W = v/ejiPA, where ejipA = ^ — vPq 
and Pq ~ GqGq is the dynamical Lindhard function. In 
a self-consistent scheme, Gq is replaced by the interact- 
ing Green function G. The Fock self-energy thus be- 
comes Yqw ~ GW, hence the name of the approxima- 
tion. Formally, the GW approximation can be obtained 
by Hubbard-Stratonovich decoupling the Coulomb inter- 
action V via an auxiliary bosonic field </> characterized by 
the propagator W ^ {4'4')- This amounts to replacing 
the electron-electron interaction by the indirect interac- 
tion of two electrons mediated by a boson described by 
(j). The first-order self-energy diagram in the expansion 
of this electron-boson interaction is E ~ GW. 

The standard derivation of the GW approximation re- 
lies on a truncation of Hedin's equations?^ where the 
three-legged vertex A = 1 -I- ^GGA is set to unity. 
In the following, we will derive the GW approximation 
for our lattice model in a diagrammatic way based on 
the four-legged vertices of standard perturbation theory, 
for both free energy functionals, that is, both choices 
of the HS decoupling. In the HS-V approach, the lat- 
tice action of Eq. ([T7| contains two interaction ver- 
tices: a local electron-electron interaction Urii^ni^ and 
a local electron-boson interaction i(f>ini. Consequently, 
the perturbation expansion of the Green's function will 
contain two types of bare interaction vertices, namely 
ro?(ri,T2,T3,r4)yfe; = USijkiSi-fSjiSk^SiiS{Ti ~T2)S{t3 - 

TA)5iT2-T3) and r|,°^(Ti, r2, ra)^^ = iSijkS{Ti-'T2)S{T3~ 
T2), which we will represent as shown in Figs. [1] and [5] 
We will suppose that we can perform the expansions sep- 
arately and then sum the two results (which is an ap- 
proximation since there could well be sequences of inter- 



actions with alternating Fc? and 



cb . 



In the HS-UV 

(0) 



approach, there is only the electron-boson vertex F, 
The vertex Fo? will lead to contributions that are not 
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FIG. 1: Diagrammatic expansion of the electron-electron ver- 
tex. From left to right, top to bottom: Bare electron-electron 
interactionvertex Fee' • Fully boldified Ecc. Expansion of the 
full electron-electron vertex, Fee. 






^AA/ 



FIG. 2: Diagrammatic expansion for an electron-boson ver- 
tex. From left to right, top to bottom: Bare electron-boson 
interaction vertex f'^'. Fully boldified Ecb. Fully boldified 
n. Expansion of the full electron-boson vertex, Feb. 



present in the HS-UV decoupling scheme. The pertur- 
bation expansion in powers of Fc? yields a series of self- 
energy diagrams, whose lines are non-interacting Green's 
functions Go. Since some higher-order diagrams contain 
"self-energy insertions" , the number of diagrams can be 
reduced by "boldifying" the lines, namely by replacing 
Go by the interacting Green's function G. Subsequently, 
the number of diagrams can be further reduced by re- 
grouping the interaction vertices into a "boldified" vertex 
Fee pictured in Fig. [T] Thus the electron-electron part of 
the self-energy (beyond the Hartree self-energy) can be 
described (Fig. [T]) by the exact expression^ 



En. - f1?gggf. 



(28) 



with bold propagators G and bold vertices Fee. Note 
that only the right vertex is boldified to avoid double- 
counting. This vertex is also called the fully reducible 
vertex. We stress that this electron-electron contribu- 
tion to the electronic self-energy is absent in the HS-UV 
approach, since in this scheme there is no longer any 
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electron-electron interaction vertex after the HS decou- 
pling. 

Similarly, the expansion of the partition function in 
powers of can be simplified by boldifying the Gq 
and Wq {= V or w"') lines and the vertices, leading to 
graphs of the form displayed in Fig. [21 corresponding to 
the expressions 



2ri°^GGTeb. 



(29) 
(30) 



The expansion for the bold electron-boson vertex Feb is 
also pictured in Fig. [51 

The usual GW approximation consists in truncating 
the expansion of the electron-boson vertex function Feb 



after its first term, namely by taking Feb ~ T 



(0) 
cb 



id. 



This simplification, which amounts to neglecting the so- 
called "vertex corrections", yields the familiar expres- 
sions E^J^ = -GW and U'^'^ = 2GG. For the HS-UV 
decoupling, S = Ecb and thus 



^H^-uv — ~GW, 



(31) 



For the HS-V decoupling, there is a second contribution 
coming from the electron-electron vertex. If one approx- 
imates Foe ~ ^cc (as in Ref. J^), one gets = 
-U^GGG, and hence 



^HS-V 



-GD - U^GGG. 



(32) 



At this point, a few remarks are in order: The two ap- 
proximations (151]) and are not equivalent. Making 
the lowest-order approximation on the electron-electron 
vertex is a stronger assumption than truncating the 
electron-boson vertex. In the HS-V approach, the series 
of diagrams corresponding to Eq. ([5^ contains the ring 
of "bubbles" made up of G and (which contains only 
the off-site repulsion V), plus a second order diagram in 
U. In contrast, the diagrams in contain the 

ring of bubbles made up of G and Vij (which contains U 
and V) . Put differently, it not only comprises the off-site 
interaction V to all orders (at the RPA level), but also 
the on-site interaction U to all orders (at the RPA level). 
We thus expect the HS-UV scheme to be better poised to 
capture non-local effects arising from V and U, while the 
HS-V scheme will probably give non-trivial contributions 
only in parameter regimes where V plays the dominant 
role. While we benchmark both approaches in the results 
section, we therefore focus in the following discussion of 
the combined GW-I-DMFT scheme on the formulation in 
terms of the Almbladh functional. A similar combina- 
tion based on the vj/y functional is possible, leading to 
a combination of the GD-l-self-consistent 2nd order per- 
turbation theory expression of Eq. p2l) with dynamical 
mean field theory, as described in Ref. 14. We refer to this 
combination in the following as "GD-t-SOPT-}-DMFT" . 



C. GW+DMFT approach 

As already hinted at before, EDMFT and GW are com- 
plementary approximate schemes: EDMFT provides a 
good description of local correlations, while GW cap- 
tures longer-range correlations and in particular long- 
range screening. Therefore, combining both approxima- 
tions appears promising. The GW-t-DMFT approachi^iii 
consists in making an approximation on 'i/[Gij,Wij] by 
decomposing it in the following way: 

* « *^^M^^[G,„ H/,,] + ^^Zoc{G^j,W.j] (33) 



where ^^Zoc 



GW 



loc • 



While EDMFT will generate the series of local self- 
energy diagrams up to infinite order, the non-local con- 
tributions to ^ will be generated in a perturbative way 
by the non-local part of GW diagrams, thus avoiding 
double-counting. 

In the limit of infinite dimensions, non-local diagrams 
vanish. Thus, the effect (if any) of the non-local con- 
tributions is expected to occur only as the dimension is 
lowered. The proximity to a phase transition should also 
enhance spatial fiuctuations. In principle, the GW con- 
tribution should nonetheless remain a correction to the 
DMFT part, which justifies why ^' nonioc [G^ , Wij] can be 
treated on a perturbative level, while ^'[Gii, Wa] is eval- 
uated exactly. 

The approximate electronic self-energy will be given by 



'{l — 6ij)Yiij where 



Sloe and = 



GW 



The 1 — Sij factor ensures the substraction of the local 
part of the GW self-energy. Analogous expressions hold 
for Hy . 

This approach is very general. In the specific case of 
the extended Hubbard model, one can expect the GW 
contribution to be significant as one approaches an in- 
stability in the charge sector, namely close to the charge- 
ordering transition. The GW diagrams can in principle 
be replaced by any other perturbative diagrammatic cor- 
rection, corresponding to decoupling the interaction in 
other channels. 



IV. NUMERICAL IMPLEMENTATION 

A. Solution of the EDMFT impurity problem 

The impurity models (|23l) and (15^ can be solved ef- 
ficiently using the hybridization-expansion continuous- 
time quantum Monte-Carlo solver (CTQMC-hyb).^^ The 
formalism has been previously derived using a Hamilto- 
nian representation of the impurity model. ^'^'^^ Here, we 
discuss an alternative derivation based on the effective 
action, focussing on the case of action (|^ . A CTQMC- 
hyb simulation consists of sampling configurations repre- 
senting specific time-sequences of "hybridization events" , 
with weight proportional to the determinant of a matrix 
of hybridization functions. The perturbation expansion 
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of the partition function Z and the summation of dia- grams with identical operator sequences leads to 
I 



^ n 

{n„}=0 cr 



X 



1 rP fP fP 

J-— drl dr'l... dr^l d^^DetA, 
[na') Jo Jo Jo Jo 

J i^[c*,c]e-^-Tnc:(rTK(rn...c:«Jc.«J, 



(34) 



where {Aa)ij — Ao-(Ti — rj) is the hybridization func- 
tion evaluated for the time- difference between annihila- 
tion operator i and creation operator j and 



Sat = j dTdT'na{T)lAaa'{T - T')n„f{T') 

+ fdTY,(^Arr{dr-fi)c,{T) (35) 



represents the interaction and chemical potential con- 
tributions of action (|24p . The interaction U can al- 
ways be split into a delta-function contribution and a 
non-singular contribution, U[t)^ ^, = Ub{T){\ — 5a, ct') + 
I'(t)ct,(t'. (In the HS-V scheme, this separation is al- 
ready explicit.) The last factor of Eq. (p4| can be easily 
evaluated since the time evolution operators are diag- 
onal in the occupation number basis. In the segment 
representation^^ each imaginary time interval with occu- 
pation no- = 1 is marked by a segment, and the last factor 
can (up to a permutation sign) be written as w^wi/Wx) 
with = efJ-ih+h) and W[j = g^^'"™''"? . Here, stands 
for the total length of segments of spin cr, while /overlap is 
the total overlap between segment of opposite spin (see 
illustration of a segment configuration in Fig. The 
retarded interaction contributes a factor 

^ g-5 E„i,„2 Jo'' '^-^i /o'' '*^2X'(ri-r2)„i,„2"<'l(-^i)"<'2(^2) 

(36) 



/overlap 




FIG. 3: (Color online). Illustration of a Monte Carlo config- 
uration in the segment representation. The top figure corre- 
sponds to a configuration with one spin-up and one spin-down 
segment, each representing a time-interval marked by a cre- 
ation operator (empty circle) and an annihilation operator 
(full circle). The overlap of the segments, corresponding to 
the width of the hashed region, yields the weight due to the in- 
stantaneous interaction U. The retarded interaction ©(r — r') 
is represented by the blue curved lines. The bottom panel rep- 
resents the weight of this configuration after integration of the 
retarded interaction over the segments. The dashed blue lines 
correspond to the interaction K{fi ~fj) between creation and 
annihilation operators. 



tegral in Eq. ([55)1 thus yields 



where {fco-} represents the collection of segments of spin 
a. 

Let us now define a function K{t) such that K"{t) — 
V{t) for < t < /3 and K{0) = K{fi) =0. if is /3- 
periodic and symmetric around t — (i/2. It has a slope 
discontinuity at zero, so that the second derivative also 
gives a delta-function contribution. Carrying out the in- 



+ if'(0)(/t + /;)+2if'(0)/overlap ■ (37) 

Using the fact that K{t) is an even function, we can 
write Inwv = J2^>J s^SJ[K{f, - t,) - if (0)] + if'(0)(/t + 
l^) + 2iir'(0) /overlap where the time- arguments of the hy- 
bridization events (creation and annihilation operators) 
are now ordered as < ri < f 2 < . . . < /3 and s is -1-1 for 
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a creation operator and —1 for an annihilation operator. 

We conclude that the retarded part of the interaction 
I?(t — r') results in a retarded "interaction" between all 
pairs of impurity creation and annihilation operators, as 
well as a shift of the instantaneous interaction = 
U — 2i^'(0) and a shift of the chemical potential /i — > 
/i = /i + iir'(O). If one writes the interaction term in 
terms of density fluctuations, ^ / / n(T)2?(r — T')n{T') 
with n = n — (n), the only change induced in the weight is 
yet another shift of the chemical potential jl — /i+i^'(0)— 
2(n)if'(0). The retarded interactions can be evaluated 
at negligible computational cost since the calculation of 
this contribution for a local update is 0{n) (where n 
is the number of operators), while the evaluation of a 
determinant ratio is OirP'). 

In practice, the local bosonic propagator W\oc ~ 
{T(j){T)4>{0)) needed in Eq. (p6)) is deduced from the 
connected charge-charge correlation function xioc — 
{Th{T)n{0)) via the relation 

Wioc =U~ UxiocU. (38) 

Indeed, using Eq. (|22|) , W\oc can be reexpressed as W\oc = 
-2|Jii#. The chain rule = -^^^^ and Eq. ^ 

give = — 5XI0C + ^7 and hence one arrives at 

Eq. ((38|) . An analogous expression holds for Dioc- 

B. Self-consistency 

The EDMFT+GW scheme is generally expected to 
work well if the non-local GW contributions to the self- 
energy is a relatively small correction to the local self- 
energy computed by EDMFT. It thus makes sense to 
first obtain a reasonable guess of the final solution by 
applying EDMFT only, and then compute the non-local 
correction and study its effect on the properties of the 
system. Following this observation, we implemented the 
EDMFT+GW scheme as follows: for a given U and V, 
we (i) obtain a converged EDMFT solution, (ii) take 
the EDMFT solution as the starting point for a self- 
consistent EDMFT-I-GW calculation, and (iii) stop when 
local and non-local observables have converged. 

This is not the only possible combination of EDMFT 
with GW, although the final result of the self-consistent 
scheme should not depend on the starting point, provided 
the scheme converges. For instance, one can choose to ini- 
tiate the scheme by computing and 11'^^ from non- 
interacting propagators Go{k, iuJn) — [«w„ + — e(A:)]~^ 
and Wo (A:, ivn) — Vk- Yet, these propagators yield a very 
large GW polarization owing to their metallic character 
(they correspond to the U = 0, V = case). Especially 
for the regions of interest here {V close to Vc and finite 
U), this large polarization is far from the expected solu- 
tion. Indeed, one can observe that when taking an insu- 
lating G and e.g. W = v as inputs for GW, the resulting 
polarization is small compared to the local polarization 

Hloc- 



C. Analytical continuation 

The non-trivial structures of the frequency-dependent 
interaction result in additional features in the local spec- 
tral function A{uj) = —^ImG[uj + irf). For example, the 
case ImP(w) — — A'^ {^{'^ ^ ^0) ~ + "^o))? studied in 
Ref. 1^1, corresponds to the Holstein-Hubbard model, for 
which the local spectral function is expected to display 
plasmonic peaks at multiples of the "plasmon" frequency 
wo--^ However, the commonly- used MaxEnt algorithm^ 
tends to smooth out high-energy features and there- 
fore a dedicated scheme must be implemented to re- 
cover the sought-after features. We proceed as proposed 
in Ref. 19: (a) From U^ivn), we compute the bosonic 

function Bir) = cxp (e„^o ^^^^^^-^ {e"^- l)) , its 

Fourier transform B{ivn) and, using a Fade procedure,— 
its spectral function B(lu); (b) From G{t), we compute an 
auxiliary function Gaux(''') — G{t) / B{t) and use MaxEnt 
to obtain Aaux(w); (c) Finally, we compute the spectral 
function as the convolution 

pea _|_ g-/3<^ 

A[u.) = J ^ _,-,./ au.(a; - e). 

(39) 



V. SUMMARY OF THE COMPUTATIONAL 
SCHEME 

The computational scheme can be summarized as fol- 
lows for the HS-UV (resp. HS-V) scheme: 

1. Start with an initial guess for S(fc, iuj) and Tl{k, iv): 
for instance, S = and 11 = (non-interacting 
limit). 

2. Lattice Green functions: Compute G{k, iw) and 
W{k,ii') (resp. D{k,iv)) via Dyson's equation 
with Vk (resp. w^') as the bare interaction. 

3. EDMFT self-consistency: Extract Gioc(»w) = 
^j,G(fc,icj) and W\oc{w) = X]fe^(^'*^) (I'esp. 
D\oc) and use Eqs. JM]) and ^ to find Q{iu) 
and U{iv) (resp. Viiv)). 

4. Impurity solver: Compute Gioc(t) and Xioc(''') 
(resp. Xioc('^))' as well as Wioc = U-Ux\ocU (resp. 
D\oc = "D — Vxi^^V) . From these, extract the self- 
energies Sloe = G^^ ~ G^^^ and IIioc =U^^ - W{^1 
(resp. (ny)ioc = V-^ - D-^l). 

5. GW-I-DMFT step (optional): 

(a) Compute 

nGW'(fc,T) = 2E,G(g,r)G((7 - k-r) 

and E«'^(fc,r) = -Y.qG{q,T)W{k - 
(7, r), where = W — v (resp. 

EGH'(fc,r) = ^E,G{q,r)D{k ~ q,T) ~ 
C^'E,G(<z,r)nGi^(g-fc,r)). 
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(b) Extract non-local parts from GW: 



and U^^i^^^^{k,iiy) 



(c) Combine Sioc(«'^) and iw) 
into iw), as well as Hiodii^) and 

nS:riocai(fc,*^) into n(fc,zz.). 

6. Go back to step 2 until convergence. 

In a pure EDMFT scheme, steps 5 (a, b, c) are skipped. 
Note that the form of the GW self-energy, GW^, is cho- 
sen such that no double counting of the Hartree energy 
occurs. Figure 2] gives an overview of the implementation 
of the GW+DMFT scheme. 




Gi„, Wioc 



g u 

IMPURITY SOLVER 
Ximp 



^imp 



Hi 



mp 




= Himp + n; 



FIG. 4: Computational scheme (HS-UV decoupling). 



VI. RESULTS 

In this section, we present numerical results for the U- 
V Hubbard model using the different approximate for- 
malisms discussed in the previous sections. We solve the 
impurity problems using the CTQMC-hyb method. Un- 
less otherwise stated, the computations are performed 
at inverse temperature /3-D — 100. The /c-sums are dis- 
cretized in the irreducible Brillouin zone on a 80 x 80 
grid, while the imaginary time correlation functions are 
measured on a grid of A'^ = 1000 equally spaced points. 
Up to 40 EDMFT steps are required to reach convergence 
close to the Mott transition. 



A. Phase diagram 

Figure [5] shows the phase diagram in the space of the 
parameters U and V for the two decoupling schemes HS- 
UV and HS-V. The top panel shows the EDMFT result, 
the bottom panel corresponds to GW-I-DMFT. There are 
three phases: (i) a Fermi liquid (FL) metal at small U 
and small F, (ii) a charge-ordered (CO) insulator at small 



U and large V , and (iii) a Mott insulating (MI) phase at 
large U and small V. 

The phase boundary to the charge-ordered phase has 
been located by approaching the phase transition from 
below Vc- The phase transition corresponds to a diverg- 
ing charge susceptibility, namely to the formation of a 
pole in the Fourier transform x(fc,aj) of Xij{^ ~ = 
d{ni{t)) /dhj{t'), where hj{t) is a probe field. Specifi- 
cally, the charge-ordering transition will be signaled by a 
divergence at Q = (tt, tt) and u = since the probe field 
for this phase is hi{t) — he'"'^^'. Using the action one 
can easily show that Xijit~'t') = (^i(*)^i(^'))- Recalling 
that W = V — vxv, we find the exact relation 



U{k,Lu) 



1 - VkTl{k,uj) 



(40) 



for the HS-UV scheme. Similarly, for the HS-V scheme, 
Xv can be computed from D — w"' — v'^''xvv"'^ or 
Xv = ^ny/(l — ti"'ny). This shows that the transi- 
tion also corresponds to the appearance of a pole in the 
fully-screened interaction W{k,ih'n), and provides a rig- 
orous definition of Vc for HS-UV and HS-V, respectively: 

1 - (;7-4V;)n[A: = (7r,7r),w = 0] = 0, (41) 
l + WcUvik^ {tt,'k),uj = 0] = 0. (42) 

On the other hand, the phase boundary between the 
metal and the Mott insulator is signaled by a vanish- 
ing spectral weight at the Fermi level, which is related 
to the imaginary-time Green's function by A{uj = 0) = 
lim^_>.oo §G(|). The curvature of the FL-MI phase 
boundary shows that increasing the nearest-neighbor re- 
pulsion V makes the system more metallic. 

Within EDMFT, both decoupling schemes yield very 
similar phase diagrams. The boundary of the charge- 
ordered phase is characterized by two main regimes: for 
U < Uc ~ 2.5, dVc/dU « 1/4, which is the prediction 
of mean-field studies. For U > Uc, dVc/dU sa 2. The 
transition between the two regimes is marked by a kink. 
This kink also coincides with the point where the charge- 
ordered critical line meets the Mott critical line Uc{V). 
The latter is only weakly dependent on V . In the temper- 
ature range [3D g [25, 100], the phase diagram does not 
depend much on temperature. The critical value Vc{U) 
for t/ < [/c is substantially larger than its naive mean- 
field estimate. We note, however, that EDMFT may 
overestimate the effect of the local interaction, leaving 
the true value of Vc at an intermediate value. 

The effect of the GW contribution to the phase dia- 
gram depends on the decoupling scheme. For HS-UV, 
GW does not have any influence of the phase bound- 
aries, while in HS-V, GW substantially lowers the FL- 
CO phase boundary. This has the following origin: the 
HS-V scheme resums the diagrammatic series for V and 
for U separately (and treats U only to second order), 
whereas the HS-UV scheme resums both terms simul- 
taneously. HS-UV is thus better poised to capture the 
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FIG. 5: Phase diagram of the U-V Hubbard model for HS-UV 
(red), and HS-V (green) decoupUng. The top panel show the 
EDMFT result, the bottom panel the GW+EDMFT result. 



competition between the localizing term U and the delo- 
calizing term V . That GW in this scheme does not alter 
the phase boundaries should therefore come as no sur- 
prise: it merely shows that the local (EDMFT) physics 
alone fixes the critical value of the non-local interaction, 
and suggests that the HS-V decoupling underestimates 
Vc- For this reason, we will henceforth restrict most of 
our attention to the HS-UV scheme. 

Fig. [6] plots the results for ImGioc(jw„), ReWioc(iw„), 
ImEioc(ii^ri) and ReHioc(«w„) corresponding to the 
EDMFT simulation at U — 2.2 and various values of 
V. As V grows, |ImG(«a;o)| increases and |Sioc(ia;o)| de- 
creases, which indicates that the system becomes more 
metallic as a result of screening by V . Indeed, the screen- 
ing effect can be quantified by the static values of the 
fully screened interaction Vl^(O) = W\oc{ivQ) and of the 
partially screened interaction C/(0) = Uiivo), which are 
plotted in Fig. [51 Nearest-neighbor repulsion V induces 
a screening of the on-site Hubbard U, which becomes 
C/(0). When V increases, 1^(0) and C/(0) get smaller 
and smaller, resulting in a more metallic behavior. The 
transition to the charge ordered insulator occurs around 
when the cost of doubloon formation vanishes, i.e. when 
W{0) = 0. 

The last panel of Fig. [B] shows the polarization 
nioc(i^'n) (which is the local bosonic self-energy). 
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FIG. 6: Imaginary-frequency data for the EDMFT calcula- 
tions at 17 = 2.2 and indicated values of V (HS-UV scheme). 
Panel (a): ImG'ioc(ii^n)- Panel (b): ImSioc(i<^n)- Panel (c): 
ReWioc{ii^n)- Panel (d): ReIlioc{ii^n) ■ 
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B. Screening 

1. Screened interaction and screening frequency 

The off-site interactions transiate into an effective re- 
tarded interaction at tire ievel of tire impurity action, as 
made apparent in Eq. (j23p . Tire frequency dependent 
local interactions in the HS-UV formalism are now de- 
scribed by U{uj) (or U + 'D{uj) in the HS-V decoupUng 
scheme). An important observable is W\oc^ which is the 
fully-screened interaction, as opposed to lA, which ex- 
cludes local screening effects. We have decided to fo- 
cus more specifically on W\oc , which we have analytically 
continued to real frequencies using a Pade scheme*S3, The 
shape of Wioc(w) a.tU — 2.2 and various V is displayed in 
Fig. [71 Wioc(w) has the typical shape of a screened inter- 
action: the real part features two distinct energy scales, a 
bare interaction Woo = W\oc{'^ oo) — U at high ener- 
gies and a screened interaction W{0) = Wioc(w = 0) < U 
at low energies, separated by a screening frequency wq. 
Its Kramers-Kronig-conjugated imaginary part has most 
of its spectral weight concentrated around wq. 





FIG. 7: Influence of V on the effective local interaction for 
U — 2.2 (HS-UV scheme). Upper panel: HeWiod^)- Lower 
panel: ImWiod^)- 

In the following, we will mainly focus on three param- 



eters to describe screening effects: (i) the value of the 
static screened interaction VF(0), (ii) the strength A of 
this screening, which we will define later, and (iii) the 
screening frequency cjq. 

W{0) is the effective fully-screened interaction between 
two electrons on the same lattice site. Its evolution across 
the U-V plane for the HS-UV scheme is illustrated in the 
middle and lower panels of Fig.[8l W{0) decreases for in- 
creasing V, and drops to zero as V approaches Vc- This 
is intuitively easy to understand: the critical line cor- 
responds to the locus where the cost Wioc(w = 0) for 
the formation of doublons vanishes. The lower panel of 
Fig. [5]) shows that the screening of the local interaction 
is much more efficient in the metallic phase than in the 
Mott insulator. While the screening is especially large in 
the metallic phase, it is weak and weakly ^-dependent 
in the Mott insulating phase (at least for low to interme- 
diate values of V). 

Let us emphasize that there are screening effects even 
when V — in the metallic phase (middle panel, red 
curve with crosses). This shows that in a EDMFT de- 
scription of the simple Hubbard model {V = 0), there 
is a screening of the static U by the local polarization 
caused by U itself, provided one uses the HS-UV decou- 
pling scheme. In the HS-V method, the screening comes 
only from V, as D originates from the HS decoupling of 
the nearest-neighbor interaction only. The upper panel of 
Fig. [8] shows U{uj = 0) for comparison. This corresponds 
to the local static interactions without polarization ef- 
fects. As expected, W(w = 0) > 1^(0) since the local 
polarization further screens the local interaction. 

Motivated by the Hamiltonian representation of the 
impurity model with dynamically screened interaction 
(see Appendix [Cj we define the strength of the screen- 
ing by the parameter A = W | dulmU . It follows 



from Eq. (jC5j) that A cx ^JJ2p where Ap is the cou- 
pling of the harmonic oscillator with frequency uip to the 
charge on the impurity, tiq. Therefore, A can be inter- 
preted as the strength of the coupling to the charge fluc- 
tuations. Its dependence on U and V is presented in 
the upper panel of Fig. [HI Except in the vicinity of the 
phase transition to the charge ordered phase, A is propor- 
tional to V. The square of the proportionality constant 
(see inset) decreases with increasing U , and exhibits two 
regimes separated by a kink at U ^ Uc- for U < Uc, 



dU 



\dvJ 



-2.2, otherwise ^ 



KdV > 



-0.34. 



This again indicates that screening in the Mott insulator 
is radically different from screening in the Fermi-liquid 
metal. 

A relevant question is in which parameter regime 
a model with a static screened interaction pro- 
vides a reasonable approximation. A useful quan- 
tity to investigate in this context is the screening fre- 
quency Wo, whose precise determination is a some- 
what tricky task owing to the Fade procedure's inac- 
curacy. Instead of measuring loq as the minimum of 
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FIG. 8: Static screened interactions in the HS-UV scheme. 
Upper panel: partially screened U{u = 0) as a function of U. 
Middle panel: fully screened VKioc(<j-' = 0) as a function of U. 
Lower panel: W^ioc(i^ = 0) as a function of V. 



ImWloc(i^), we have found it more meaningful to de- 
termine it from the first moment of ImWioc(w): ujq ~ 
(cj) = dujujlmWioci^)/ dujlmWioci^), whose de- 
pendence on U and V in the HS-UV scheme is pre- 
sented in Fig. |9] (middle panel). This figure shows that 
the screening frequency is only weakly dependent on the 
nearest-neighbor interaction V. On the other hand, the 
larger the bare interaction U — Uao , the larger the screen- 



ing frequency. The lower panel shows that the screening 
frequency (for V = 0) has two distinct regimes separated 
by a kink sX U = Uc- 



2. Screening frequency from linearized DMFT 

The U -dependence of the screening frequency may be 
traced back to the form of the one-particle local spec- 
trum. As can be seen from Eq. ([55]) . for V = the 
frequency-dependence of W^ioc (and thus the value of wq) 
is inherited from the charge-charge correlation function 
Xioc- An analytical estimate for the poles of this func- 
tion can be calculated by means of a simple approxi- 
mation named linearized DMFT.^"* In this method, the 
impurity problem is approximated by the coupling of the 
correlated impurity to a single uncorrelated bath level 
describing the hybridization of the impurity to the lat- 
tice degrees of freedom. This simplified version of the 
impurity problem allows for the explicit calculation of 
the local Green's function)^ 



i=i ^ 



W -f 



(43) 



with 



1 = \ {^U^ + MV^,^T ^fu^Tm^ , (44) 

,(45) 



as well as 02 — ^ — ai. The hybridization strength's 

dependence on U is given by Vhyb = i^/i-^/l — W^/U^ 
(see Ref. [U for details). Uc denotes the critical U for 
the Mott transition. In the estimate below, we will use 
the value computed within EDMFT, Uc = 2.5. 

In the absence of vertex corrections, the corresponding 
charge-charge correlations function can be computed as 
Xioc = —26*1006100, leading to the expression 



Xioc(w) = -2 



2ale2 



4aia2e3 



(46) 



The six (3 x 2) poles are defined as ei — 2ei, €2 = 2e2 
and 63 — El + 62- 

These poles, displayed in the lower panel of Fig. [HI 
correspond to the transitions allowed in the various cor- 
relation regimes, namely: in the low-correlation limit, 
only transitions within the quasi-particle peak are possi- 
ble. As correlations increase, the appearance of Hubbard 
bands enable additional transitions from the lower Hub- 
bard band to the unoccupied states of the quasi-particle 
peak, and from the occupied states of the quasi-particle 
peak to the upper Hubbard band. In the strong corre- 
lation regime, finally, the only possible transitions are 
those between the lower and the upper Hubbard band. 
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FIG. 9: Screening parameters. Upper panel: A as a function 
of V (HS-UV scheme). Inset: {j^f as a function of U. Mid- 
dle panel: (uj) as a function of V. Lower panel: i'CD-XioiF^'^'^ 
in the {pJ,U) plane for V = Q. The poles have been artifi- 
cially broadened by an imaginary factor rj = 0.01. Black line: 
EDMFT result for {to) as a function of U (for V = 0). 
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FIG. 10: Local spectral function (HS-UV scheme). Upper 
panel: U = 2.2. Lower panel: U = 3.5. 



V tend to smooth out the incoherent Hubbard bands 
and transfer some spectral weight to the quasi-particle 
peak and into the gap region between the quasi-particle 
peak and the Hubbard bands (Fig. [TUl upper panel). 
This behavior is consistent with the imaginary-frequency 
data showing a more metallic behavior as V increases. 
The effects are more dramatic in the strongly correlated 
regime {U > Uc), where one observes new features in the 
local spectral function, as shown in the lower panel of 
Fig. [TOl In addition to the two Hubbard bands located 
a,t oj — ±U/2, the spectral function has two symmet- 
ric satellites at lu — zt{U/2 + ujq), whose spectral weight 
grows with V. The position of the peaks comes from 
the convolution in Eq. (1391) . since Aaux contains spectral 
weight at ±U/2 and B contains weight at ±ujq. 



C. Spectral Properties 

The screening effects coming from U and V have some 
impact on the local spectral function. In the weakly- 
correlated regime {U < Uc), the non-local interactions 



D. Momentum dependence 



We now turn to the analysis of the GW+EDMFT re- 
sults. 
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1. Non-local self-energies 

Figure E] displays iwo) and 

nnoniocai(^' ^'^o) in the metallic phase near the charge 
ordering transition. These quantities vanish in the 
limit of large dimensions and are thus neglected in the 
EDMFT treatment. The GW contribution to the imag- 
inary part of the electronic self-energy S is negligible 
with respect to the local self-energy (for instance, at 
U = 2 and V = 0.6, ImSioc(«wo) = ~0.18, compared to 
a non-local GW self-energy < 0.001). This holds across 
the Fermi-liquid phase and the Mott insulating phase. 
The real part of Sj^J^j^j.^^ is relatively large away from 
the EDMFT Fermi surface, but does not alter this Fermi 
surface. On the other hand, the non-local polarization is 
comparable to its local counterpart (nioc(wo) = —0.39 
iovU = 2,V = 0.6). 
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FIG. 11: GW Non-local contribution in the HS-UV scheme 
at ?7 = 2 and V = .4 (left) or 1/ = .64 (right). Upper panel: 
ImS^^;(,c(fc, ia;o). Middle panel: KeT,^^i^^{k,iujo). Lower 
panel: Ren^Zioc{k,ii^o)- 

This should have an impact on the phase diagram if 
one recalls the criterion of Eq. (|42|) . However, the local 
observables are also modified in the self-consistent calcu- 
lation, as will be described in the next subsection, which 
prevents a direct prediction of the effect of the non-local 
terms. As shown in Fig. [SI the effect of GW depends on 
the decoupling scheme. For the HS-UV scheme, the GW 
contribution has a negligible influence on the phase dia- 



gram. For the HS-V scheme, the GW contribution has 
a large effect on the phase boundary between the metal- 
lic phase and the charge-ordered phase. The non-local 
polarization, peaked at fc = (tt, tt), enhances nesting ef- 
fects and leads to a substantially lower Vc compared to 
EDMFT. 



2. Influence of the self-consistency on local observables 

As already hinted at, the self-consistency condition 
leads to an adjustment of the local quantities to the non- 
local self-energies. This is illustrated in Fig. [T^l which 
displays the local spectral function after EDMFT con- 
vergence and after GW-I-EDMFT convergence, and in 
Fig. [131 which illustrates the convergence from EDMFT 
to EDMFT+GW. 




-2-101-2-101 



00 (0 

FIG. 12: Effects of GW on the local spectral function ^ioc(i^) 
(HS-UV scheme, U = 2). 

The EDMFT+GW local spectral function features 
more pronounced Hubbard bands than the EDMFT spec- 
trum. This is a sign of increased local correlations. In- 
spection of the imaginary frequency data corroborates 
this observation: Wiodi^ — 0) is enhanced with respect 
to the EDMFT result, indicating that the local interac- 
tions are stronger. Likewise, |ImGioc(j'j-'o)| is reduced by 
GW. 

These observations mean that the local quantities have 
become more "insulating" in character as a result of 
the addition of the non-local GW self-energy. This 
can be interpreted in the following way: contrary to 
the EDMFT case, where all the screening and correla- 
tion effects are absorbed into the local self-energy, in 
EDMFT-I-GW some of these effects are now carried by 
the non-local components. Specifically, GW non-local 
self-energies carry important screening effects owing to 
the very nature of the GW approximation. This leads to 
a redistribution of the screening between local and non- 
local observables: local observables become less screened, 
and thus more correlated. 
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FIG. 13: Convergence properties for !7 = 2, V = 0.4 (HS- 
UV scheme). EDMFT corresponds to the converged EDMFT 
resuh, EDMFT+GWi denotes the observables corresponding 
the zth iteration with non-local self-energies, with EDMFT 
as starting input. Upper panel: ImGioc(i<^n). Lower panel: 
ReVKioc(«i'n). 



3. Convergence properties and self-consistency 

Figure [TH shows the converged ImGioc and ReVFioc 
for the three self-consistent schemes: EDMFT alone, 
EDMFT-I-GW and GW alone (with the EDMFT-^GW 
result as a starting point). As expected, GW is the 
most metallic in character and the corresponding spec- 
trum does not have Hubbard bands. Interestingly, 
GW-HEDMFT is not a kind of "average" between GW 
and EDMFT. It exhibits stronger correlation effects than 
both GW and EDMFT. 




FIG. 14: Influence of the self-consistency fov U = 2, V = 0.4 
(HS-UV scheme). Upper panel: ImG'ioc(iWn). Middle panel: 
IieWioc{ivn) ■ Lower panel: ylioc(t*j) (obtained from MaxEnt 
continuation) . 



VII. CONCLUSIONS 

We have presented an application of the EDMFT and 
GW-hDMFT methods to the single-band extended Hub- 
bard model. In a first step, the two formalisms have been 
reviewed in detail: We have discussed the construction 
of the free energy functional, and compared two differ- 
ent flavors of such functionals that have been proposed 



in the literature, corresponding to two distinct decou- 
plings of the interaction term, and leading respectively to 
the GW+DMFT and GD-f SOPT+DMFT approaches. 
We have presented the details of our implementation 
of a fully self-consistent GW-I-DMFT scheme based on 
a numerically exact continuous-time quantum Monte- 
Carlo solver adapted for frequency-dependent local inter- 
actions. The investigation of the frequency-dependence 
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of these interactions for parameters ranging from weak 
to strong coupling shows that the J7- dependence of the 
local screening frequency reflects the form of the local 
one-particle spectrum. We have investigated the spec- 
tral properties of the extended Hubbard model within 
three self-consistent schemes, namely EDMFT, GW, and 
GW-I-DMFT. The nearest-neighbor repulsion V leads, in 
the Mott insulator, to high-energy satellites in the local 
spectra. 

The GW-I-DMFT calculations demonstrate that the 
non-local contributions to the self-energy coming from 
the GW diagrams are quite small in the case of the ex- 
tended Hubbard model. In view of the strong momen- 
tum dependence observed in self-energies obtained from 
cluster DMFT calculations for the two-dimensional Hub- 
bard model as one approaches the Mott transitionj^iS 
our results confirm the importance of spin fiuctuations, 
suggesting that further non-local diagrams have to be 
considered in order to capture the dominant fiuctuations 
in the extended Hubbard model. 

The model calculations presented in this paper can 
be straightforwardly extended to the multi-orbital case 
and to additional, longer-range matrix elements of the 
screened interaction, paving the way for realistic first- 
principles material calculations. It is worthwhile to note 
that the GW-I-DMFT method and its variations are a 
computationally cheap way of incorporating the leading 
vertex contribution - in the form of the EDMFT self- 
energy - into the description of a solid, and to introduce 
some spatial fluctuations through a perturbative scheme. 
This contrasts with methods involving an explicit com- 
putation of the vertex functionsii^ whose implementation 
for simple model systems is already a formidable chal- 
lenge. 

In real materials, further degrees of freedom, stem- 
ming for example from the multi-band nature and lig- 
and orbitals, lead to a renormalisation and/or frequency- 
dependence of the parameters in the low-energy descrip- 
tion. Relatively weak, but non-local correlation effects 
are expected to play a dominant role in the case of ex- 
tended ligand or higher-lying empty states, thus pro- 
viding an additional motivation for a combination of 
GW and DMFT. Indeed, GW provides an accurate and 
comparatively inexpensive description of the screening 
from "imcorrelated" bands, making the application of the 
GW-I-DMFT method to electronic structure calculations 
for realistic solids highly promising. 
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Appendix A: Extended Hubbard models from the 
First Principles 

The Hamiltonian H of the extended Hubbard model 
is supposed to provide the low-energy physics of cor- 
related materials. H can be regarded as an effective 
Hamiltonian resulting from a "downfolding procedure" , 
which is expressed in some localized Wannier basis. The 
downfolding procedure is akin to a renormalization group 
transformation which, starting from all the bands result- 
ing from a LDA calculation, produces an effective model 
for the bands in a low-energy window by integrating out 
the remaining bands. In this process, the bare Coulomb 
interaction v{r,r') = 47re^/|r — r'| is transformed into 
to a frequency-dependent partially screened interaction 
Wr{r,r' which acts in the low-energy subspace. 
In principle, Wr is computed as Wr — v{l — vPr)~^, 
where Pr is the polarization computed when transitions 
within the effective model have been removed. The ma- 
trix elements of this interaction in the Wannier basis are 



y,,fezH = J d\d\'^*{r)4,*{r')Wr{r,r',uj)Mr)Mr'), 

(Al) 

where denotes a Wannier orbital centered at site i. 

Model ([ij involves three approximations on the 
above matrix elements: (i) the frequency-dependence of 
V^jkiiuj) is neglected: Vtjki = Vijki{uJ = 0), (ii) the in- 
teraction is restricted to density-density terms Vijki ~ 
VijijSikSji = Vij, and (iii) only the on-site matrix ele- 
ment U = Vii and the nearest-neighbor matrix element 
V = Vij (with i and j nearest neighbors) are retained. 
The last assumption is valid only if Wr{r — r',uj = 0) de- 
cays rapidly in space^^^ The neglected non-site-diagonal 
parts of the electron-electron interactions such as, for 
instance, the bond-charge-bond-charge matrix elements 
W = Vijji, are believed to be small in usual solids-^ 



Appendix B: Derivation of the single-site EDMFT 
action using the cavity method 

In the following we use the cavity method^ to derive 
the EDMFT action, Eq. (HH), and the EDMFT self- 
consistency equations which flx Q and U. To this end, 
let us focus on a given site (denoted by the index 0) 
and split the lattice action (Eq. ((T7)) ) into three parts: 
S = So + -I- A^* where 5*0 denotes the action of the 
site 0, S^°^ the action of the lattice with site removed 
(the lattice with a "cavity" at site 0) and AS the remain- 
ing part: 
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So 



AS 



5(0) 



+ icpouo + -(f>o[v \oo(po}, 



J ^'^ I X! '^0'^ ~ ^) ^O"- 



(Bl) 
(B2) 

(B3) 



Defining rjia- = UqCq^ and j.^ = [v ^]io</'o, we can write / T)^o[c*,Ci,(j)i]( 



-(So + S'^' + AS) 



/Z. We can express 



AS = IodT[~j:^^,^^irj*c,, + clr,,)+j:^^,j,^,y the ^action as 5eff = 5o - f^K* ,j7,^j,] +^const, where 
such that 7]ia and ji can be regarded as sources ^iVi iVidi] = In / 2?5^o[Cj , c^, (?!)i]e '■^ is the gen- 

of correlation functions for the effective action erating functional of connected correlation functions of 



of the site 0, defined by e-'5'=ffK.'=«'*«l/^o 



the cavity)^ 



g: 



(0) 



Sjn (n) . . . 6j,„ {Tn)6jj, (t{) . . . Sjj^ «) ' 



(B4) 
(B5) 



An explicit expression for fl is thus 

oo „ 

m,VuM = E E J dri.-. dr'^inl in) . ■ . vn (n . . . <) 

GO „ 



(B6) 



n=l ii...i„ j'i...j„ 



The DMFT approximation consists in approximating 
by its infinite-dimensional limit. In this limit, the hop- 
ping t between sites must be scaled as t/y/z (with z — 2d) 
in order to keep a finite kinetic energy, while V must be 
scaled as V/z in order to keep the Hartree energy cor- 
responding to the nearest-neighbor interaction finite4i 
As a consequence of taking this limit, all terms of 
order n > 1 in Eq. ((B6)) vanish, so that Q,^'^^^^ = 

/ dTdT'clir) (- U^t.^Gf (r - r')) co(r') -|- 

/ drdr'Mr) (E., vi,\j,'w';^^ (r - r')) Mr')- We 
thus arrive at the DMFT effective action of Eq. ([2T|) 
(where we have dropped the index in order to simplify 



I 

the notation) with 



-YUotjoGf^ii^n), (B7) 



^1 Y^^.-l^-lM/W 



00 



^m^^W^r('^n). (B8) 



Furthermore, in the limit of infinite dimensions, the 
cavity Green's function is related to the lattice Green's 
function through G^f = Gij - GmGoj/Goo and W^j^ = 
Wij — WioWoj /Wqo, which is shown by considering the 
paths contributing to Gij (Wij) and not to g|^^ {W^^^) 
(see Ref. i for more details). This allows us to write 
^iotjoGf'\ after Fourier-transforming, as 



E4Gfc(jw«) 



ekGkiiuJn) 



(B9) 
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At this point a second approximation is made: the 
self-energies are assumed to be fc-independent, namely: 
io;) w I]ioc(jw) and n(fc,w) nioc(w)- This also 
becomes exact in the c? — >■ cxd limit. ®^ As a consequence, 
we can define the densities of states p{e) — 5{e — e^) 



and /9'(e) 
(|B9]) as 



which allows us to rewrite 



dep{e)e 



(BIO) 



where C = iw„ + /i — I]ioc(i'^n)- The same expression 
holds for the screened interaction, with p p' and C, — > 
C = [w^^loo - nioc(w„). 

Using the following identities for Hilbert transforms: 



J —i 



dep{e)f_ 
= dep{t)€ 



c 



dep{e)e 



= -i + c 



(Bll) 



(B12) 



we obtain the self-consistency relations (|25p and ([26| . 

Equations (PT|) . (^51) and form a closed set of equa- 
tions: S'^*^^^, once solved, yields Sioc and IIioc, which 
gives updated G and which can in turn be used to 
solve the effective local problem again until convergence 
is reached. 



Appendix C: Hamiltonian formulation of the 
impurity problem 



Some properties of action (I24p are more easily under- 
stood in terms of its hamiltonian representation. The 
first two lines correspond to an Anderson impurity model 



Ha 



IM 



E 



p 

(CI) 



On the other hand, the retarded effective 



'i—'P iu),n—e 

interaction can be generated by coupling the impurity to 
a bath of bosonic modes described by the Hamiltonian 



^^^(</.2+n2)+^Apno0p, (C2) 



with (j)p = ^ (bp + bl) and Xlp = {bp - 6],) . 

Using the identity 11^ = —{dr(t)p{T))'^/ujp, this can be 
written in an action formulation as 



Jboson = — ^.^VpyiVrn) I ;^7^ ) (Pp{~lVm) 



+ Xp(pp{iiy,n)n{ 



2uj. 



p 



(C3) 



Integrating out the bosonic degrees of freedom leads to 



^boson = "q X] ^'^(^^^'^) < A: 

m y. P 



2 '2ujp 



Defining V[iu^) = J ^lmP(c^) u.Jv^-u- ^i^h 



n(-Wm). 

(C4) 



Iml?(a;) = — tt 2, ^pH'^ ~ ^p) 



p 



(C5) 



and Fourier-transforming 'D{ivn) yields the retarded in- 
teraction in Eq. (IM)) . 

The retarded interaction may thus be regarded as 
stemming from the coupling to a bath of harmonic os- 
cillators labeled by the index p, with frequency Wp and 
coupling strength Ap (as already emphasized in Ref. Tl'). 
The effective interaction mediated by these auxiliary de- 
grees of freedom is proportional to the squared coupling 



describing an impurity (c, ct) coupled to a bath of non- strength times the free-phonon Green's functio: 
interacting fermionic levels (a^, a|,, energy £p ). Here, Dl^iv^) = -{TMr)MO)) 
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Ha = c]yCa and n — -I- n^. The connection be- 
tween Eqs. (|24l) and (jCip is given by Q~'^{iu)m) = i^jJm + 
p — A(za;m) and the hybridization function /S.{iu)„i) = 



T- — v5 — T- Note that in 
complete analogy to the fermionic hybridization function 
A(aj), the frequency dependent interaction I'(aj) is de- 
termined self-consistently. 
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